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We study particle current in a recently proposed model for coherent quantum transport. In 
this model a system connected to mesoscopic Fermi reservoirs (meso-reservoir) is driven out of 
equilibrium by the action of super reservoirs thermalized to prescribed temperatures and chemical 
potentials by a simple dissipative mechanism described by the Lindblad equation. We compare 
exact (numerical) results with theoretical expectations based on the Landauer formula. 
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I. INTRODUCTION 



Particle current trough a coherent mesoscopic conductor connected at its left and right hand side to reservoirs 
is usually described in the non-interacting case by a formula, due to Landauer, based on the following physical 
picture: Electrons in the left (right) reservoir which are Fermi distributed with chemical potential (M-r) ^^nd 
inverse temperature /3l {Pr) can come close to the conductor and feed a scattering state that can transmit it to the 
right (left) reservoir. All possible dissipative processes such as thermalization occur in the reservoirs while the system 
formed by the conductor and the leads is assumed to be coherent. The probability of being transmitted is a property 
of the conductor connected to the leads which is treated as a scattering system. In this picture, the probability that 
an outgoing electron comes back to the conductor before being thermalized is neglected, the contact is said to be 
reflectionless. This description of the non-equilibrium steady state (NESS) current through a finite system has been 
rigorously proved in some particular limiting situations [T], such as infinite reservoirs, but several difficulties prevail 
the understanding of non-equilibrium states in general and the description of current in more general situations, for 
instance in the case of interacting particles. Two frameworks are usually considered to study these open quantum 
systems: one deals with the properties of the state of the total (infinite) system [H [5] where reservoirs are explicitly 
considered as part of the system. The other is based on the master equation of the reduced density operator, obtained 
by tracing out the reservoirs' degrees of freedom and is better suited to be applied to different systems and to compute 
explicitly some averaged NESS properties, at the price of several approximations such as e.g., Born-Markov (see e.g. 

IS)- . 

In this paper we explore particle current in a model where we mimic the leads that connect the reservoirs with the 
system, as a finite non-interacting system with a finite number of levels (which we call meso-reservoir). The reser- 
voirs (called here super-reservoirs) are modeled by local Lindblad operators representing the effect that a Markovian 
macroscopic reservoirs have over the meso-reservoirs. In Sec.|Tl]we introduce the model and briefly review the method 
we use to solve it. In Sec. |III| we analyze the particle current operator and indicate the quantities that should be 
computed for a full description of the current. In Sec. Ill A we briefly present the Landauer formula that is expected 
to apply in some appropriate limits to our model and in Sec. |IIIB| we analyze the numerical results we obtained with 
our model and compare them with the current predicted by the Landauer formula, validating the applicability of our 



model but also going beyond by computing the full probability distribution function (PDF) of the current. In Sec. IV 
we present some conclusions and discuss interesting perspectives of our study. 



II. DESCRIPTION OF THE MODEL 



We consider a one-dimensional quantum chain of spinless fermions coupled at its boundaries to meso-reservoirs 
comprising a flnitc number of spinless fermions with wave number k (k Cz {1, . . . ,K}). The Hamiltonian of the total 
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system can be written as H = Hs + + Hr + V , where 

N-l N 

is the Hamiltonian of the chain with {tj^ the nearest neighbor hopping, U the onsite potential and Cj,c^ the annihi- 
lation/creation operator for the spinless fermions on the site j of the chain (conductor). The chain interacts through 
the term 

K 

^ = E {^k^lL^^i + ^'f + (h.c), (2) 

k=l 

with the meso-reservoirs Ha — ^fe^Ia^fca- Here a = {L, i?} denotes the left and right meso-reservoir. They 

share the same spectrum with a constant density of states 6^ in the band [-Emin, -^max] described by eu = Og^k — fcg) 
and afe,Q,a^ ^ are annihilation/creation operator of the left and right meso-reservoirs. The system is coupled to the 
leads only at the extreme sites of the chain with coupling strength that we choose fc- independent 114) = Va- 

We assume that the density matrix of the chain - meso-reservoirs system evolves according to the many-body 
Lindblad equation 



d 
dt 



k. 



where m G {1,2} and Lk,a,i = ^ Fa{ek))aka, ifc,Q,2 = VlPai^kjata operators representing the coupling 

of the meso-reservoirs to the super-reservoirs, Fa{e) = (^e^aie-fj^a) ^_ i^-i gj-^ Permi distributions, with inverse temper- 
atures Pa and chemical potentials pa, and [•, •] and {•, •} denote the commutator and anti-commutator, respectively. 
The parameter 7 determines the strength of the coupling to the super-reservoirs and to keep the model as simple as 
possible we take it constant. The form of the Lindblad dissipators is such that in the absence of coupling to the chain 
(i.e. Va =0), when the meso-reservoir is only coupled to the super- reservoir, the former is in an equilibrium state 
described by Fermi distribution. [71 [S] . 

To analyze our model we use the formalism developed in [9| . There it is shown that the spectrum of the evolution 
superoperator is given in terms of the eigenvalues Sj (so-called rapidities) of a matrix X which in our case is given by 



1 7 

X = —-Hi^ay + —\ONxK OnxN 



xK 0, 
OkxK Orxn 




where O^xj and Ej denote i x j zero matrix and j x j unit matrix, ay is the Pauli matrix, and is a 
matrix which defines the quadratic form of the Hamiltonian, a,s H = SHd in terms of fermionic operators 
= {aiL, ■ ■ ■ aKL,Ci, • • • , CN,aiji, ■ ■ ■ aKu}- 
The NESS average of a quadratic observable like djdi is given 9 in terms of the solution Z of the Lyapunov equation 
X'^Z + ZX = Mi with M, = -|diag{TOiL, • • ■ ,mKL;^\xN-, mui, ■ ■ ■ .rriKR} ® cTy and uika = 'y{'2Fa{sk) - 1} as 
follows: Consider the change of variables W2j-i = dj + dj, W2j = i{dj — d^), the NESS average of the quadratic 
observable WjWk is determined by the matrix Z through the relation (wjWk) — Sj^k ~ ^i'^j.k- Wick's theorem can 
be used to obtain expectations of higher-order observables and in fact, the full probability distribution for these 
expectation values in some cases. 

III. PARTICLE CURRENT 

The operator representing the current flowing from the fc-th level of the meso-reservoir to the chain is given by 

^^^'fclalci-clafe), (4) 

while the current trough the site I of the chain is 

Ji = iti{c]ci+i - cJ+iQ)- (5) 



3 



In the steady state the average current is conserved in this model (HIS] and thus (J/) is independent of I. Moreover, 
if we define the current from the left meso-reservoir as J = Y^^=i3k^ have that (J;) = (J). 
It is not difficult to note that current satisfies 

Jr-{'Ll% (6) 
I Ji if "-even ^ ' 

with J° 1. Now we are in a position to compute the full non equilibrium current distribution in terms of (J;) and 
{Jf). For this wc consider the generating function 



^ {iky 



n=0 



(7) 



which, using Eq.(|6]), gives 



»fcJ,\^^_W^^V:^PQg^^^A:iiAjsinfci. (8) 



if J tf t, 



The probability distribution p{Ji) is the inverse Fourier transform of {e thus we get 



which is normalized. 

We note that normality and positivity of probability lead an interesting inequality: < < 1 — An 
equivalent result holds for the current from the fc-level of the meso-reservoir to the system: 



These expressions are expected because ±ti and are the possible eigenvalues of the operator J; (similarly for 
jj^), but they show that (J;) and (Jf) contains all the information about the current. We will study these quantities 
numerically, thus we need to solve the above mentioned Lyapunov equation and note that in the Wj variables (we use 
the primed variables to indicate indices as they appear in variable d, i.e. if j is a site in the chain, then j' = K + j) 

Jl = i^-^{w2l'~lW2l' + l + W2VW2l'+2), (11) 

thus 



{J]) — '2tj{Z2j'-l,2j' + l + Z2j',2j'+2) (12) 

and 

jf = ^(1 + W2j'-lW2j'W2j' + lW2j'+2)- (13) 

Wick's theorem implies 

(Jj) = ^(l - 16(Z2j'-l,2/+2Z2j',2j' + l - Z2/-l,2i' + lZ2j',2j'+2 + Z2/-l,2j' Z2/ + l,2j'+2)) • (14) 

To simplify the discussion, in what follows we reduce the number of parameters by assuming constant hopping and 
onsite energy ti = t and Ui = U . Moreover we fix 7 such that ^ > Va- In that case, we showed |2j that the transport 
quantities become roughly independent of 7. In fact, in that case, the coupling super-reservoir - meso-reservoir is 
stronger than that of the system (chain) - meso-reservoir. Then the meso-reservoir are driven to a near equilibrium 
state weakly dependent of 7. We explore now the behavior of the current as a function of flp, Wq, t and contrast the 
observation with expectations based on the Landauer formula. Qualitative explanation of the current behavior is also 
provided. 
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A. The Landauer Formula 



The Landauer formula |10j provides an almost explicit expression for the NESS average current as a function of the 
parameters of the system. In units where e = 1 and h = I it reads 



(J) 



1 

271 



where T(a;) = tv\r i^[lj)G'^ {(jj)T if {uj)G {lj)] is the transmission probability written here in terms of 

G±(a;) 



(15) 



(16) 



the retarded and advanced Green function of the system connected to the leads and of —Tat/^ the imaginary part of 
the self-energy E+. 

The self energies have only terms at the boundaries of the chains JD] i.e., (E^)„m = cr^(5„m^n,b, where 6 = 1 if 
a = L and b — N ii a — R and 



K 



fc=i 



(17) 



Recalling that both leads have the same spectrum, we assume a constant density of lead states 1/Oq in the range 
[£'min , ^'max] , thus Ta{Lo) — 2Trv'^/0Q IS independent of uj inside the interval and zero otherwise. For the real part of 
the self-energy we have the principal value integral 



A„(a.) = ^P 
zvr 



UJ - 



(18) 



The Landauer formula is expected to hold when the leads have a dense and wide spectrum. Therefore, we restrict 
ourselves to the case that -Emin ^ ~t and t <C -Bmax, the so called wide-band limit, where Aq(w) can be neglected. 
The transmission coefficient is then T{uj) = TlTjjIG^j^Iu})]'^ and we need to compute the wide-band limit retarded 
Green function 



/uj + i^f -1 



G+{u;)- 







a; -1 
-1 ■■. 



V 







(19) 



-1 u + i^ J 



Note that in the previous expression we have set [7 = which sets the energy axis origin and < = 1 which sets the 
energy scale. Thanks to a recursion relation, this matrix can be inverted jllj and one explicitly finds the relevant 
element of the Green function 



I I 

{2uJ+^^{^,+Tn))EiJ ^("1) 



k^N-2-2k(N-2-k 



k^N-4-2k /'JV-4-fc\ 



(20) 



where [xj is the largest integer smaller than x. In the next section we compute numerically the integral in Eq.(15) 
and compare with the results obtained in our model. 



B. Numerical results 



In Fig. [T]we depict in blue and red two Fermi distributions with K = 50 levels and parameters ^^=4, /3i = 3 and 
fiR — —4, /3fl = 3 respectively. In the middle (brown) the spectrum of a chain with 7 sites and t = 3. The width of the 
chain energy band is 4i with 7 levels inside and is centered around U — 0. From this picture, we expect that decreasing 
the width A/i of the populated energy interval [fiR,^^L] is equivalent to increase the width of the conduction band t. 
This is confirmed in the left panel of Fig. [2] where we show that the current is roughly independent of t for 2t < fi^ 
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FIG. 1: Schematic plot of the left and right Fermi distributions and the levels of the middle chain with N = 7 sites. 
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FIG. 2; Current versus hopping f for 7 = 0.1, U = 0, Pl ^ = 5 , vl = vr ^ 0.05, K = 200 and TV = 10. In the left panel 
(black) fiL = 4, hr — —4 and (red) /il ~ 10, hr = —10. Right panel: Same parameters except; t — 1, vl = vr = 0.03. 



and decreases with t for 2t > (black dots), when the conduction band extends beyond the region populated by 
electrons in the reservoirs. The red dots are obtained for a larger A/i for which the conduction band is always inside 
the populated region. In the rest of our numerical examples we set [/ = and t ^ 1. Analogously, in the right panel 
of Fig. [2] we show that for fixed t the current grows linearly with A/i and saturates at A/i > 4t. 

What is perhaps more interesting is the scaling of the current with F = ^L^a/i^L + ^n)- In Fig. |3]we plot the 
current (J) as a function of the coupling to the right lead, showing that the main trend of the current is (J) = F. In 
the inset we show that there are deviations to this law. 

We have analyzed this behavior using the Landauer formula, that allows deeper analytical exploration. Since 
temperature is very low we take it exactly zero in both reservoirs, thus the Landauer formula is (J) — 
^l^r/'2'T^ /y^j" tr[G+(a;)G~(aj)]. Changing the representation of the Green matrix from position basis to the en- 
ergy basis (of the isolated chain) one can see that it is possible to approximate (isolated resonances approxima- 
tion) the integrand as a superposition of N Lorentzians each located at 2i cos[j7r/(A^ -I- 1)] and with a width 
(Tl + Fij) sin^[j7r/(iV + 1)]/(A^ -I- 1), j = 1, . . .iV along the energy axis. If i < /i^, the integral can be extended 
from —00 to 00 because the Green function decays exponentially outside the energy band of the chain, thus the result 
(J) = F is obtained. 

Now we can also explore fluctuations and compute (J^). In Figj4^, we computed for the same parameters than in 
Figj2^ the quantity 2( J^)/<^ — 1. We see that as soon as the current decreases because some modes of the chain goes 
out of the populated energy region, the fluctuation increases. As a function of other parameters like v^, or N (data 
not shown), the average (J^) does not have important changes (see Fig jib). 



IV. CONCLUSIONS 



We showed that in the wide-band limit, the numerical results found in our model indeed correspond to what is 
expected on the basis of the Landauer formula, a formula which is usually interpreted as if the reservoirs were always 
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<J/r) 




FIG. 3: Black points are obtained by our model with 7 = 0.1, Pl = Pr = 5, = —4, = 4:-, N = 10, U = 0, t = 1, vl = 
0.03, K = 200 and -Emin = —20, Smax = 20. The red dots are obtained by a numerical computation of the Landauer formula 



with T — 0, and the continuous line is (J) = T — 2n 



The inset is (J) /V. 
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FIG. 4: Panel a: 2{J^) /t^ — 1 as function of t for ^r = —4, /^l = 4 (black dots) and ^r = —10, /il = 10 (red dots) as in Fig 
III. Panel b: 2{ J^) /t^ - 1 fo r f = 3 as a function of vl = vr. In both figures the other parameters are: 7 = 0.1, Pl = Pr - 
= IQ,U = Q, K = 200 and £min = -20, £max = 20 with vl = 0.03 in panel a and f = 3, = -4, = 4 in panel b. 



in an equilibrium (gran canonical) distribution, not perturbed by the presence of the system. The Landauer formula 
emphasizes the role of the Fermi distributions of the reservoirs and provides an accurate description of the current if 
the assumption of reflectionless contacts is justified. In this respect, a very interesting relation was found in and 
proved in [3] between the current and the occupation of the meso-reservoir (J) = Y^^=i 27fe[('^fe) — Phi^k)]- This 
is an exact relation that in the appropriate limit should converge to the Landauer formula. Note that it implies 
that the occupation difference with respect to the Fermi distribution is ©(F/^qt)- It is a very interesting relation 
because it links the current, which is the fingerprint of the non equilibrium state, to the difference in distribution 
to the equilibrium case. Something similar has been found in classical systems where the fractal nature of the non 
equilibrium state is determined by the current |13| . Moreover in [2] we analyzed how Onsager reciprocity relation 
is broken in the system and found that \Lup/Lpu — 1| grows with 7 implying that despite the almost 7 independent 
value of the current, the dissipative mechanisms in the super-reservoir play an important role. A deeper study of 
these effects, that are beyond the Landauer picture, can be studied in the context of the model presented here and 
deserve further investigation. 
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